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Abstract. - We show that an electronic phase transition described by the Cahn-Hilliard equation 
has important applications to cuprate superconductors. The simulations of the local charge density 
and free energy reveal two main features: i) The segregation process creates tiny isolated regions 
with potential wells where the holes can be bound in single-particle levels, ii) The clustering 
process also gives rise to an effective two-body pairing interactions and superconducting amplitudes 
A sc (r) at low temperatures. The resulting system resembles a granular superconductor with the 
resistivity transition driven by Josephson coupling among these nanoscale grains. This approach 
reproduces the well known critical temperature transition T c (p) as function of the doping level p. 
Furthermore, the local density of states with spatial dependent gaps A(r) is due to the intragrain 
single-particle bound states that remain above T c , which characterizes the pseudogap phase and 
reproduces many measurements. 
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Introduction. — The origin of the superconducting 
gap associated with the superconducting state and the re- 
lation to the pseudogap above the transition temperature 
(T c ) remains one of the central questions in high-T c re- 
search [TH3]. It is a matter of debate whether this prob- 
lem is connected with the charge inhomogeneities observed 
in many experiments [4--6. and with the position depen- 
dent energy gap A(r) measured by Scanning Tunneling 
Microscopy (STM) [THIS]. In this letter we suggest that 
these two problems are strongly related and their solution 
must be worked out together. 

There are steadily accumulating evidences that the 
charge distribution is microscopically inhomogeneous in 
the CuO-i planes of high temperature superconductors 
(HTSC). The most well known charge inhomogeneity is 
the charge stripe structure (HIS]. Nuclear quadruple reso- 
nance (NQR) have measured two signals corresponding to 
two (low and high) densities around the average doping [B] , 
that increases as the temperature is decreased. This evolu- 
tion with temperature can be an evidence of an electronic 
phase separation (EPS). More recently, NQR experiments 
indicated that the charge inhomogeneity in the planes is 
correlated with the O dopant atoms in YBa^Cu^Oy sys- 
tem [14] . Also a granular fractal microstructure of intersti- 
tial O was recently measured in LaiCu^y [15j Similarly, 



recent STM data have revealed one of the most puzzling 
property of cuprates, the non-uniform energy gaps A(f) 
on the length scale of nanometers [7hT5]. The derived 
LDOS have generally two forms: one type with smaller 
gaps and well defined peaks and other with larger gaps and 
ill-defined peaks [5] . Furthermore, some of these gaps re- 
main well above the superconducting critical temperature 
T c (p) [91H0]. These two types of gaps [5] may occur due to 
the existence of two energy scales that have been observed 
on electronic Raman scattering measurements [TB], STM 
data [TThT5] . Angle Resolved Photon Emission (ARPES) 
[17] and combined STM- ARPES [18]. However, the origin 
and even the existence of this two-gap picture is still a 
matter of debate [TW^D] . 

In order to obtain an unified interpretation to all of 
these experiments we study an EPS transition that gen- 
erates regions of low and high densities. Such transition 
may be driven by the lower free energy of undoped anti- 
ferromagnetic (AF) regions [21] (intrinsic) or by the out 
of plane dopants [T3] (extrinsic origin) like oxygen inter- 
stitials [T5]. As p increases, the Coulomb repulsion among 
the charges in the local high doping regions increases the 
energy cost of the phase separation, which ceases the tran- 
sition in the overdoped region, in agreement with the dis- 
appearance of the local AF fluctuations [35]. This clus- 
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tering phenomenon gives rise to a potential that bunch 
up the holes and separates low and high density phases in 
the form of small domains or grains in the Cu02 planes. 
This EPS process can be regarded as due to an effective 
attractive potential, which is used in the self-consistent 
calculations to obtain the intragrain superconducting am- 
plitudes. In this approach the resistivity transition tem- 
perature T c (p) occurs when the Josephson energy Ej(p) 
among the separated regions is equal to KbT c (p). 




Fig. 1: (color online) The density map simulation of the inho- 
mogeneous charge density on a grid with 105 x 105 sites. The 
white square indicates where the superconducting calculations 
are made. For the p = 0.16 compound, the local densities pi 
along the line with a black trace are shown on the top of the 
figure. 



The Potential from the Charge Inhomogeneities. 

— Phase separation is a very general phenomenon in 
which a structurally and chemically homogeneous system 
shows instability toward a disordered composition [23.. As 
discussed above, it has been recognized to be very impor- 
tant to cuprate superconductors [53] but there is not a 
theory to describe the degree of disorder as function of 
doping and temperature. Consequently, to deal with this 
problem, we use the general theory of Cahn-Hilliard (CH) 
[25 27 adapted to the cuprates. The starting point is 
the EPS temperature Tps(p), taken close to the upper 
pseudogap that produces an anomaly measured in some 
experiments [TH3]. The transition order parameter is the 
difference between the local and the average charge den- 
sity u(p,i,T) ee (p(i,T) -p)/p. Clearly u(p,$,T) = 
corresponds to the homogeneous system above Tpsip). 
Then the typical Ginzburg-Landau free energy functional 
in terms of u(p, i, T) is given by 



f(u) = -e 2 \Vu\ 2 + V(u,T). 



(1) 




Fig. 2: (color online) The simulation of the potential 
V(p, i, T). The values on 42 sites at the white line are shown in 
the top inset to demonstrate the potential wells with average 
barrier ~ V g b where single-particle bound states are formed. 
The low inset shows schematically an example of the V (p, i, T) 
minimization with the holes (black dots) clustering by the CH 
diffusion as an effective two-body hole attraction. 



size of the boundaries between the low and high density 
phases 26,27 . The CH equation can be written [23J in the 
form of a continuity equation of the local density of free 
energy /, d t u = —V.J, with the current J = A/V(<5/ /5u), 
where M is the mobility or the charge transport coefficient 
that sets the time scale. Therefore, 



du 



2„,3\ 



-MV 2 (e 2 \7 2 u - A 2 {T)u + B 2 u 



(2) 



Where the potential V(u,T) = -A 2 {T)u 2 /2+£V /4+..., 
A 2 (T) — a(Tps(p)—T), a and B are constants, e gives the 



In Fig.©, we show a typical simulation of the density 
map with the two (hole-rich and hole-poor) phases given 
by different colors. 

The potential V(p,i, T) isolates the hole-rich and hole- 
poor regions forming grain boundaries between these two 
phases [57]. V(p,i,T) reaches its minimum value at the 
low and high equilibrium densities as it is demonstrated 
by the simulations shown in Fig.© and Fig.©. Our most 
important finding is that V (p, i, T) produces two different 
effects: first, it creates potential wells as it is shown in 
the top inset by the values of V(p,i, T) on 42 sites along 
the white straight line. The potential barriers among the 
grains can be approximated by V g b(p,T) and confines the 
holes, generating local or intragrain single-particle bound 
states. The signature of these single-particle bound states 
will appear in the local density of states together with 
the superconducting gap from the two-body attraction. 
Second, and more important, to minimize the local free 
energy, the holes move to equilibrium positions in a similar 
fashion as if they attract themselves. This is schematically 
illustrated in the low inset of Fig.© where a) represents 
a homogeneous system with p = 0.25 (one hole at each 
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four sites) and b) the motion toward clusters formation 
of low (pi = 0) and high densities (pj = 0.50), according 
to the CH diffusion process. Clearly, this movement of 
holes can be regarded as originated from an effective two- 
body attraction, like the spin-spin exchange interaction, 
that arises from the Pauli principle and the electrons wave 
functions. We do not know the exact strength of this 
potential but it should scale with the potential barrier 
V g b(p,T) between the two phases. We discuss below that 
this effective potential can give rise to superconducting 
pairs. 

The Self-consistent Calculations. The calcu- 
lated u(i,T) (or p(xi,T)) on square lattices like the 
42 x 42 shown in Fig.fTT]) is used as the initial input 
and it is maintained fixed throughout the self-consistent 
Bogoliubov-deGennes (BdG) calculations. The pairing 
potential V g b(p, T) is adjusted to match the average local 
density of states (LDOS) measured by low temperature 
STM on 0.11 < p < 0.19 Bi2212 compounds [SJ, namely, 
V gb {p,T = 0) = (0.7 - 2.36 x p) (eV). Starting with an 
extended Hubbard Hamiltonian, the BdG equations are 



k a \ / u „( Xi ) 
,A* -K*J Un(xi), 



'un(xi) 

v u„(xi) 



(3) 



These equations, defined in detail in Refs. [27H3T] . are 
solved self-consistently. u n , v n and E n > are respec- 
tively the eigenvectors and eigenvalues. As mentioned, 
the hole clustering process is mainly along the Cu — O 
bonds, and with the Coulomb repulsion, it favours the d- 
wave pairing. Thus, the d-wave pairing amplitudes are 
given by 



A d (xi 



5) + w*(x. i )M„(x i 



+5)] tanh 



E,, 



2k B T' 

and the inhomogeneous hole density is given by 



(4) 



As the top inset of Fig.© shows, the different local 
densities regions are bounded by the potential barriers » 
V g b- As the temperature goes down, the barriers increases 
forming an metal-insulator-metal junction. Consequently, 
the electronic structure in the CuO — 2 planes becomes 
similar to that of a granular superconductor [35]. In this 
way the superconducting transition occurs in two steps: 
first by the appearing of intra-grain superconductivity in 
the grains and than by Josephson coupling with phase 
locking among all the hole-rich and hole-poor regions at a 
lower temperature. 
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T|«n(x i )| 2 / n + l«n(x i )| 2 (l- /„)], (5) 
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Fig. 3: Top panel, the average A<j(p), and in the inset, the 
schematic single particle and superconducting energy levels 
present in each nanoscopic grain. In the low panel, the thermal 
energy ksT and the Josephson coupling among superconduct- 
ing grains Ej(p,T) for some selected doping values as function 
of T. The intersections give T c (p), as plotted in the inset. 



where /„ is the Fermi function. We stop the self-consistent 
calculations only when all p(xj) converges to the CH den- 
sity map shown in Fig. ([I]). 

Due to the Cu d-orbitals and the strong on site Coulomb 
repulsion, we calculate only the intra-grain d-wave su- 
perconducting gap Ad(i,T). The effect of the tempera- 
ture in the potential is taken into account by V g b(T) ~ 
(I - (T/Tps) 1 - 5 , as demonstrated by CH [25] near T PS , 
and consequently, Ad(i,T) — > at a single temperature 
T*(p). The fact that all gaps vanishes at the same tem- 
perature although the values of A d (i,T — 0) can be very 
different is a consequence of the self-consistent mean field 
approach. 



By using the theory of granular superconductors [33J 
to these electronic grains and the calculated average su- 
perconducting amplitudes Af(T,p) = £)f A d (T,i,p)/N, 
where N is the total number of sites, we can estimate the 
values of T c (p) . 



Ej(p,T) 



nhA a /(T lP ) 
4,e 2 R n 



tanh{ 



(T,P) 



2K B T C 



)■ 



(6) 



Where A^" (T, p) is the average of the local superconduct- 
ing gaps A d (i,p,T) on a N X N (N = 28, 36 and 42) 
square lattice which is plotted in the top panel of Fig.([3]). 
The inset shows the schematically single particle levels at 
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a shallow puddle whose walls are proportional to V g b- R n 
is the normal resistance of a given compound, which is 
proportional to the planar resistivity p ab measurements 
|34j on the La2- p Sr p Cu02 series. 

In the low panel of Fig- (J3j) , the Josephson coupling 
Ej(p, T) is plotted together with the thermal energy kpT 
whose intersection yields the critical temperature T c (p), as 
shown in the inset. The fact that T c (p) has the well known 
dome shape is due to the different behavior of A^ v {T,p) 
that decreases with p (Fig.©) and 1/R n that increases 
with p. The values are in reasonable agreement with the 
Bi2212 T c (p), as expected, since V g b was chosen to match 
the Bi2212 low temperature LDOS gaps 8 . Since, in gen- 
eral T*{p)) > T c (j>), this approach provides an interpre- 
tation to the superconducting amplitude and the measured 
quasiparticles dispersion in the normal phase |19U20| . 

The LDOS and superconducting gap. In the 

BdG approach, the symmetric local density of states 
(LDOS) is proportional to the spectral function [35J and 
may be written as 

Ni(T,V gb ,eV) = ^[K(x 4 )| 2 + Mx,)| 2 ] x 

n 

[f' n (eV-E n ) + f' n (eV + E n )}.(7) 

The prime is the derivative with respect to the argu- 
ment. u n ,v n and E n are respectively the eigenvectors 
and eigenvalues of the BdG matrix equation [271125] , 
/„ is the Fermi function and V is the applied volt- 
age. Ni(T, Vgb, eV) = LDOS(V g b) is proportional to the 
tunneling conductance dl/dV, and we probe the effects 
of the inhomogeneous potential by examining the ratio 
LDOS{V gb / 0)/ LDOS {V gb = 0). This LDOS ratio yields 
well-defined peaks and converges to the unity at large bias, 
similar LDOS ratios calculated from STM measurements 

These features are illustrated in Fig. ([4]) for a represen- 
tative overdoped p = 0.20 compound near a grain bound- 
ary at a representative average doping hole point. At 
T = 40K, Vgb = 0.240eV and for large applied poten- 
tial difference {V > O.leV) both LDOS, LDOS{V gb + 0) 
and LDOS{V g b = 0) converge to the same values. An im- 
portant result of our calculations is that the LDOS yield 
peaks at larger bias than the values of the superconduct- 
ing amplitude Ad{i,p, T). In Fig. ((4]) we can see the LDOS 
gap Ape — 24meV (Eq.flTJ) and the much smaller su- 
perconducting gap Ad = 7.2meV (Eq.©) that is almost 
undiscernable and therefore it is marked by arrows. We 
attribute Ape to the single-particle bound states in the 
grains. In our calculations we can see that hole-poor pro- 
duces larger and not well defined peaks than those at hole- 
rich regions as noticed by McElroy et al [8j. The physical 
explanation is that, if the number of holes in a grain is 
large (hole-rich grain) they occupy higher single-particle 
levels in the local potential wells and can be more easily 
removed by the STM tip. 
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Fig. 4: The LDOS on an overdoped p = 0.20 sample for V g t = 
0.240eV and for V g i = and their ratio. The arrows show 
the values of the superconducting gap Ad — ±7.2meV , much 
smaller than the LDOS gap (~ 24meV), which is normally 
measured by ARPES and STM experiments. 



To study the pseudogap phase we perform calculations 
on an underdoped compound with p = 0.11 and T c w 
65K according Fig.©. In Fig.© we show the results at 
two representative locations in the hole-rich and hole-poor 
phases. The hole-poor region with = 0.021 has a very 
large LDOS of Ape ~ 80K and A^ = 12meV, as shown 
in the inset. At a hole-rich location {p{i) « 0.23) the 
LDOS gap is given by Apg ~ 50meV with well-defined 
peaks at low temperature, and the superconducting gap is 
Ad{T = 0) = 33meV. Perhaps due to the mean-field type 
calculations, all the gaps vanish near T s» 147K, what 
can be assigned the pseudogap temperature T*(p = 0.11). 
Likewise Fig.© the low temperature superconducting gap 
Ad produces small anomalies that are marked by arrows, 
as it was already noticed by our previous work [21] and 
detected experimentally by Kato et al [TT1I13] . 

In this scenario the pesudogap phase is formed with 
the intragrain superconducting amplitude and the single- 
particle bound states. The superconducting phase also 
contains both the local single-particle gap Apc{i,p,T) 
and the local superconducting Ad{i,p,T) but has phase 
coherent through the Josephson coupling among the 
grains. At the overdoped region T*{p) approaches T c {p) 
while their difference increases in the underdoped region. 
These two curves are shown in the derived phase diagram 
of Fig.©. We have also plotted the EPS line T PS {p) 
which opens all the phase separation process and it is close 
the observed anomaly called upper-pseudogap Timusk and 
Statt pp. 

Conclusions. — The main idea of this letter is 
the two-body electronic attractive potential and granu- 
lar structure derived from the EPS transition. We show 
that an effective hole-hole attraction appears from the 
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Fig. 6: The calculated phase diagram of cuprate superconduc- 
tors as derived from the EPS transition Tps(p) and the forma- 
tion of tiny grains. The onset of single particle bound states 
and superconducting amplitude formation occurs at T*(p). 
The dome shape T c (p) curve is due to Josephson coupling 
among the grains. These three lines are in agreement with 
the experimental results. 



Fig. 5: The p = 0.11 LDOS with V gb = 0.44eV. Top panel, 
A PG (T) at a hole-poor puddle (p< = 0.021). A PG (T = 40if) « 
80meV. The superconducting gaps A<j(T = 40K) = 12meV are 
marked by arrows. Below, the set of LDOS curves at a hole- 
rich grain, (p< w 0.20) with A PG (T = 40K) « 50meV and 
with A d (T = 0) = 33meV. In the inset A d (i,T) x T for each 
case. Both Apa(i) and A<j(i) vanish near T* = 147K. 



Ginzburg-Landau free energy and the Cahn-Hilliard simu- 
lations which reveals also a microscopic granular behavior 
with single-particle bound states (as seen in Fig. (2)). With 
the values of V g b(p) that reproduced the low temperature 
average LDOS values of McElroy et al [8] and the values of 
Tps{p) from the crossover or upper pseudogap line PUS], 
we obtain the superconducting and pseudogap phases for 
all doping p. In summary: 

i) The phase separation transition described by the 
Cahn-Hilliard equation minimizes the local free energy 
generating separated regions of different local densities. 
it) Ad and Ape have completely different nature but 
they have the same origin, namely, the inhomogeneous 
potential V g b- They are both present in the pseudogap 
and superconducting phases but distinct strength. Their 
differences were verified in the presence of strong applied 
magnetic fields and also by their distinct temperature 
behaviour in tunneling experiments [36] . 
Hi) The LDOS are dominated by the intragrain single- 
particle bound states Apa and the superconducting gap 
Ad produces only a small anomaly at low temperature 
and low bias as recently observed by some STM data 

jgiDHEig. 

iv) The observed difference in the LDOS shape, called 



"coherent" and "zero temperature pseudogap" [5], is due 
to the differences in the local single-particle bound levels 
at hole-rich and hole-poor locations respectively. 
v) The potential barriers V g h prevents phase coherence 
and, as in a granular superconductor, the resistivity 
transition occurs due to Josephson coupling among the 
intragrain superconducting regions (EqjSJ). The T c (p) 
curve is a consequence of the different behavior of 1/ R n (p) 
which increases, and A d M (p) which decreases with p. 

We gratefully acknowledge partial financial aid from 
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